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We obtained a new representation of a solution of the heat conduction 
equation with boundary condition of the third kind for a layer. The result 
is presented as a superposition of fundamental solutions for an unbounded 
system with variable coefficients, the explicit form of which is given. We 
consider the well-known problem of the evolution of the temperature field 
initially uniformly distributed in a layer. The distribution of the temperature 
field is represented in terms of the obtained functions. 
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I. INTRODUCTION 



Theoretical description of the distribution of heat in bounded solids is based on the 
solution of the heat conduction equation with the corresponding initial and boundary con- 
ditions. In the case where the heat transfer between a heated body occupying the region 
< z < L, —oo < x, y < oo and the external medium with the zero temperature proceeds 
by the Newton law, the problem is reduced to the solution the heat conduction equation 
0-H for t > 

- KATftt) = fftt) (1) 

with the boundary conditions of the third kind 

&T{r,t) 



dz 

dT(f,t) 
dz 



\ x T(f,t) = 



+ \ 2 T(r, t) = 



(2) 

z=0 



(3) 



and the initial condition 



T(f,0)=T (f) (4) 

Here, A = d 2 /dx 2 + d 2 /dy 2 + d 2 /dz 2 is the Laplace operator, f = (x,y,z), k is the thermal 
diffusivity, Ai and A2 are the relative heat transfer coefficients between the body and the 
external medium at z = and z = L, respectively, and f(r,t) is the density of thermal 
sources normalized to cp, where c and p are, respectively, the specific heat capacity and 
the density of the substance. 

The methods for solving this boundary-value problem are well-known and widely de- 
scribed in the literature both for the theory of heat conduction and for equations of 
mathematical physics 0-0]. It is worth noting that a solution of this problem is presented 
as a series in which the summation is carried out over an infinite set of eigenvalues a , where 
a are positive roots of the transcendental equation 



tanaL = a— ^ — — ^— , (5) 
cr — AiA 2 

which, in the general case, is solved numerically. In particular cases where the temperature 
of the boundary surfaces is equal to zero ( A = oo , the first boundary- value problem) or the 
heat flow through the boundaries is absent (A = 0, the second boundary- value problem), 
the roots of Eq. (|5|) are known. In these cases, the solution of the heat conduction equation is 
represented in terms of normal distributions created by instantaneous heat sources situated 
at the corresponding points of the space. This enables one, without recourse to numerical 
calculations, to carry out analytical estimations of solutions of the heat conduction equation 
for the first and second boundary-value problems for various special cases (thick layer, near 
the boundary surfaces and at the center of the layer, time temperature asymptotics, etc.). 
For arbitrary values of A, Eq. (||) was numerically investigated in detail and the results are 
presented in the form of tables and plots (see, e.g., ||). 

As far as we know, for any A > , in the literature, there is no representation of a solution 
of problem (|I])-(fD in terms of fundamental solutions for an unbounded system similar to 
those for the first and second boundary-value problems for which it is not necessary to find 
eigenvalues a of the transcendental equation (^). The possibility of this representation for 
arbitrary a is shown in the second section of the paper for Green's function. Unlike the 
known solutions of boundary- value problems of the first and second kinds, in this case, the 
fundamental solutions are in the sum with the corresponding weights that are functions of 
coordinates and time. 

As an example, in the third section, we use the obtained Green function for studying the 
evolution of an initially uniform temperature field in a layer. 

II. GREEN'S FUNCTION OF THE THIRD BOUNDARY- VALUE PROBLEM 

We consider the homogeneous heat conduction equation (|l|) with boundary conditions 
(0)-(0) for Ai = A 2 = A and the initial condition 
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T(f,0) =T (r) = 5(f-f r ), (6) 

which corresponds to a point instantaneous heat source situated at the point r ' = (x', y', z'), 
where < z' < L, and find a solution continuous in the region < z < L, — oo < x, y < oo 
for t > 0. Note that in the general case of different Ai and A2 , the expression for Green's 
function for the third boundary- value problem (|])-(@), (§]) expressed in terms of fundamental 
solutions for an unbounded system was obtained in where the Brownian motion of 
particles in a layer was investigated for various values of adsorption coefficients of particles 
by the boundary surfaces. We also noted there that the specific case Ai = A 2 should be 
considered in its own right because in going from the general expressions presented in H to 
the case Ai = A2, we have to sum expressions each term of which is divergent. 

By using methods of the operational calculus, we reduce Eq. (|l]) with initial condition 
(§) to an ordinary differential equation of the second order, the solution of which T(f,t) = 
G(f,f',t) (Green's function of the heat conduction equation) can be written as follows: 

+ exp (-*V + ihZ + (n)) {j^j + I , (7) 

where R± = r j_ — r' ± , f± = (x,y,0), r'j_ = (x',y',0), Z±(n) = Z± + 2nL, where n = 
0,±1,±2, Z± = z±z'. 

Then we use the integral representation 

exp(-k 2 z Kt) = - — —-^ / d£ exp -— + ik£ , (8) 



— 00 

change the order of integration, and take integrals over k z and £. We present the final 
result in the following form: 



Ant 



G(f,r',t)= -y E {Pt ) (Z4n),X)ew 



Z 2 (n) 



Pi+\Z + (n),\)e W (-?M)\, (9) 



where 

2|n|-l 

Pt\z, A) = 1 - (1 - 5 n0 ) £ C^(3 k (\z\,X), n = 0, ±1, ±2, . . . , 

fc=0 

PW(^, A) = 1 - £ CS&ftC*, A), n > 0, 

fc=0 
2|n|-2 

PW(^, A) = 1 - £ Cjff-i&M' A), n < 0, (10) 

/c=0 

is the binomial coefficient, S n0 is the Kronecker symbol, A = A(ztt) 1 / 2 , 

fo{z, A) = 7r 1 / 2 (2A) fe+1 F fc + ^ , fc = 0, 1, 2, ... , (11) 

F k (x) = y ^j: (exp(x 2 ) erfc x) , k = 0, 1, 2, ... , (12) 

oo 

2 /" 

erfc x = J dy exp(-y 2 ). (13) 

X 

By the direct substitution of solution (9) into Eqs. (l)-(3), we can make sure that this 
representation of Green's function satisfies both the homogeneous equation (1) and boundary 
conditions (2) and (3) (for details, see Appendix). To carry out calculations one should use 
the following recurrence relationship 

F k {x) = ^ [F k _ 2 (x) + xF fe -i(x)] , k>2 (14) 

for the function F k (x). To obtain relation (|i~4"D, we use the known recurrence relationships 
for multiply probability integrals || and representation (12) for the function F k (x). As the 
subscript k of the function F k (x) increases, the function F k (x) decreases and, for all values 
of the argument x, remains either positive (for even k and k = ) or negative (for odd 
k). The behavior of the first six functions F k (x) is shown in Fig. 1 [k — 0, 1,2,3,4,5, the 
number near a curve corresponds to the subscript k of the function F k (x)]. 
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F k (x) 
1.0 




In the limiting case A — > oo, with regard for properties of binomial coefficients, it follows 
from Eq. ( [TOP that all quantities Pj~^ and P( + > are the same and constant 

limPW(Z ± (n),A) = ±l (15) 

A— >oo 

and Eq. (|[) is reduced to the well-known solution of the first boundary-value problem for 
Green's function in a layer |i1JlO|] 



As A — > 0, [A(k£)/|Z±| <C 1], in view of Eqs. (10) and (11), we have 

ft(|Zt(n)|,A)«< l2(Kt)/j (17) 

o(A fe ), fc = l,2,... 



and 

Pt ] (Z4n), A) « 1 - 4|n|A7r 1 / 2 F {^^) , n = 0, ±1, ±2, . . . 

P«(Z + H, A) a 1 - 2(2n + l)AV/ 2 Fo (^0^) , n > 0, 

pW(Z + (n),A) « l-2(2|n| - IJAV/Vq (^p) - « < 0, (18) 



whence, for A = 



p( ± )(z ± H,o) = i 



it follows the well-known solution of the second 



boundary-value problem [|I1,|T0l 



^>--^£.{-(-^) + -(-^)) < 19 > 

In the limiting case of an infinitely thick layer (L — > oo) , we have 

^(Z_(n),A) = U 

Pl + \Z + (n), A) = 5 n0 [1 - (3 (Z+, A)] (20) 

and expression @ is reduced to the known solution of the third boundary-value problem 
for a half space ( z > 0, — oo < x, y < oo ) |3||7|.[T(| 

/ , ex p(^) r / z 2 \ / zi\ 

G(r, r , t) = \ ./ < exp -— + exp 



- 2AV/ 2 exp(A 2 + \Z + ) erfc f^pp + ^ | . (21) 



III. EVOLUTION OF A UNIFORM TEMPERATURE DISTRIBUTION 

The obtained expression ([3J) for Green's function can be used for the solution of var- 
ious boundary-value problems for a layer. As an example, we consider the simplest well- 
known problem on the evolution of the temperature in an initially uniformly heated layer 
[T (f) = To] with boundary conditions (2) and (3) for Ai = A 2 = A. The required distribu- 
tion T(f,t) is defined by the integral 



7 



oo oo oo 



T(r, t) = T J dx' J dy' Jdz' G(r, f ', t). (22) 

— oo — oo 

By virtue of Eqs. ©-(13), relation ( p2j) can be represented in the following form: 

T(f, t) = T(z, t, A) = T(^, t, A = oo) + 5T(z, t, A), (23) 

where T(z, t, A = oo) is the known solution of the first boundary-value problem 

T(z, t, A = oo) = T I £ (-l) n [erfc(nZ - + erfc ((n + 1)/ + £)] - l} (24) 

and the additional term ST(z, t, A) caused by nonzero value of the quantity A" 1 has the 
form 

{oo 2n-l 
E E Ct +1 [i? fc ,„(z) + i? fc ,„(L - z)\ 
n=l fc=0 

oo 2n ^ 

+ E E Cffii [^»(-«) + R *,n{z -L)]\. (25) 

n=0 fc=0 ) 

Here, Z = L/(4/«t) 1 / 2 , ^ = ^/(4/tt) 1 / 2 , an d 

= r — E (-2)^™^ [exp (- ((2n + Z)Z - C) 2 ) /^ m ((2n + 1)L - A) 

- exp (-(2nZ - 2 ) /5m(|2nL - z\, A)] . (26) 

It is worth to note that expressions (23)-(25) can be written at once by using Eq. (72) in 
H for the density of Brownian particles in a layer with absorbing boundaries and carrying 
out changes in the corresponding notation in it. 

As is seen from Eqs. (23)-(26), distribution of the temperature in the layer remains 
symmetric about the middle of the layer at any time t > : T(L — z,t, A) = T(z,t, A), 
which is quite natural due to the symmetry of the system under investigation about the 
plane z = L/2 and uniformity of the initial distribution of temperature. 

By collecting terms with the same indices m in Eq. (25), we can represent the quantity 
ST(z, t, A) as follows: 

6T(z, t, A) = °- V (z, A) + E [Vn(z, A) + V n (z, A)] , (27) 
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where 



2n-l 



V n (z, A) = lrn,nV m , n (z, A), 77 = 1, 2, . . . , 

m=l 
2ra 

Kifo A) = lm,nV m ,n( Z i A )' 77 = 0, 1, 2, ... , 

m=0 

K^, A) = exp (-(ei _) ) 2 ) *U4 _) , A) - exp (-(Ci +) ) 2 ) »*»(4 + \ A) 

+ exp (-(£< +) ) 2 ) ^ m (4 +) , a) - exp (-(S;\) 2 ) ^(4; } i, A), 

V^ B («, A) = exp Vm(z { +\ A) - exp [-(^if) i^O&i, A) 

+ exp (-(m) 2 ) »7m(4+n A ) - ex P (-(£i) 2 ) »7m(4+i> A), 
Vn (z, A) = (-l) n 2~ (n+1) /3 n (|z|, A), n = 0, 1, 2, . . . , (28) 

z&=2nL±z, z& = (2n-l)L±z, = 4±>/(4«t) 1/2 , ^ = 4 ±} /(4^) 1/2 . 

The quantities i m>n and / mjn are defined by the recurrence relations 

lm,n = l m -i, n + (-l)"^" 1 " 1 ^, n = 2, 3, ... ; m = 1, 2, . . . , 2n - 1, 
= -2n, 

^m,n = lm—l,n ~l~ ( 1) 2 ^271+1) 77 = 1, 2, ... J 771 = 1, 2, . . . , 271, 

io,n = 1, h,n = -2n, (29) 
which are convenient for m < n. For m > n, the following relations are more convenient: 

^m,n = lm+l,n + ( — 2) m C2™ 1 , 777 < 277 — 1, 

i <->2n— 1 

Hn-\,n — i 

lm,n = lm+l,n H (^^C^+l") 777 < 277 — 1, 

fen.n = 2 2 ". (30) 

Relation (27) is valid for any A > 0. In particular case of great values of A, which 
corresponds holding the first correction terms caused by a finite value of the relative heat 
transfer coefficient of the layer with the environment, relation (27) is simplified to the form 
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^1/21 ^ ^rn.vi 4- kivin F v v s; 





7 (n)|)(l + 2| 7 (n)|) fc 


7(n)(l + 


7 (n)|) 





a(n)|)(l + 2 


a{n)\) m 


a( 


n)(l + 


a(n)\) 



exp(-((2n+l)/-0 2 ) , (31) 



where ar(n) = ((2ra + 1)/ - OA. T(n) = (2nZ - OA, n = 0, ±1, ±2, . . . , 

2n, n > 0, f 2n - 1, n > 0, 

I -(2n + l), n < 0, -2n, n < 0. 

For a thick layer ( Z 3> 1 ), relations (23), (24), and (31) for distribution of the temperature 
can be simplified. In this case, it is sufficient to retain only terms with n = 0, 1 in sum (24) 
and the term with n = in sum (31). This yields 

T(z, t, A) « T sb (z, t, A) + T s6 (L - z, t, A) - T , (32) 

where T sb (z,t,X) and T sb (L — z,t,X) are the known distributions of the temperature in 
semibounded bodies that occupy the z > and z < L regions, respectively. In the case 
A > 1, they have the form 

exp(-£ 2 ) 



T sb (z,t,X)^T 
where 



erf£ + 



(33) 



2 * 



2 r 

erf £ = 1 — erfc£ = —q^ J dx exp(— x 2 ). 
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FIG. 2. Space and time distribution of the temperature in the layer: Bi = oo (a), 100 (b), 10 
(c), 1 (d), and 0.1 (e). 
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The space and time distribution of the normalized temperature T(z, t, Bi)/T in the layer 
calculated by the general formulae (23)-(24), (27)-(30) are displayed in Fig. 2 for various 
Biot numbers Bi = LA/2. With regard for the symmetry of the system at hand about the 
z = L/2 plane, the distribution of temperature is shown only for the < z < L/2. It 
is quite natural that the layer cools down with maximum rate in the case corresponding 
to the first-boundary problem ( Bi = oo ) when the boundaries of the layer are kept at 
zero temperature for the entire time (Fig. 2a). Note that a noninstanteneous drop in the 
temperature shown in Fig. 2a from the initial value To to zero at the initial time ( t — ) as 
the observation point approaches the boundary of the layer is not the result of calculations 
but appears due to the used uniformity of the space and time grid for the temperature field. 
For finite great values of the Biot number (Fig. 2b), the distribution of temperature inside 
the layer, in fact, does not change, whereas it takes place the smooth but fast cooling of 
the boundary of the layer. For the time of order of O.lr, where r = L 2 /Ak, the surface 
temperature is as small as 6% of its initial value. As the Biot number decreases (Bi = 10, 
Fig. 2c), the cooling of the surface in noticeably slows down and the surface temperature 
tends to zero even for t ~ r. In this case, the finite value of A influences the distribution 
of temperature deep into the layer, e.g., T(T/2,r, Bi = 10)/T(L/2, r, Bi = oo) w 1.52. For 
intermediate values of the Biot number (Fig. 2d), heat removal from the surface does not 
lead to a radical redistribution of the temperature in the central and peripheral regions of the 
layer as for Bi > 1. For example, for Bi = 1 and t — r, the ratio T(0, r, Bi)/T(T/2, r, Bi) 
is only 0.65. For small Biot numbers (Bi = 0.1, Fig. 2d), the temperature equalization 
inside the layer proceeds much faster than heat removal from the surface. This case is 
characterized by a slightly uniform distribution of temperature across the section of the 
layer [T(0, r, Bi)/T(L/2, r, Bi) 0.95] and a slow decrease in the temperature with time 
[T(L/2,r, Bi)/T wO.92]. 
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IV. CONCLUSIONS 



In the present paper, we gave a new representation for a solution of the third boundary- 
value problem of the heat conduction equation for a layer which needs no calculations of 
eigenvalues on the basis of the transcendental equation. This solution is represented in 
the form of a superposition of the fundamental solutions for an unbounded system caused 
by instantaneous heat sources situated at the points of location of sources-images of the 
layer. It is the form of representation of the solution that is widely used in the literature for 
representation of solutions of the first and second boundary- value problems. 

Since new representation (9) enables one to carry out analytic investigations not using 
graphic and table data, it can be useful for various physical problems for which, in parallel 
with numerical calculations, analytic estimations of solutions must be performed. In par- 
ticular, we can find the first additional terms caused by the account of the heat transfer of 
the layer with the environment to the known results determined under the assumption of a 
given temperature of the surface or the absence of the thermal flow through it. In the paper, 
we illustrate the possibility of finding correction terms in powers of A -1 , where A > 1, for 
a special problem of the evolution of the initially uniform temperature field in the layer. 

In addition to analytic investigations, Green's function (9) can be also useful for the 
numerical analysis of various boundary-value problems of the third order. Using recurrence 
relations (14), we significantly simplify this analysis by reducing it to the calculation of 
probability integrals and the majorant estimation of the remainder terms of the series given 
in the Appendix. 
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APPENDIX: 



We prove that solution (9) satisfies the homogeneous heat conduction equation (1) with 
boundary conditions (2) and (3). For this purpose, first, we prove the uniform convergence 
of the series in solution (9) at < z < L and t > 0. In view of the following representation 
of the quantity (3 k (\Z±(n)\, A) : 

(zL X i \ fc+1 

rar^ j • < A1 > 

which is valid for \n\ > (nt) 1 / 2 /L, we obtain that for finite values of A, for all |n| > A, 
where A > Xnt/L, the coefficients PW(Z ± (n), A) = P^^X) can be represented as follows: 

" - ir) ■ (A2 » 

As |n| — > oo, the following estimate for the coefficients holds: 

Inn pW(A)=exp(-^) <1. (A3) 

|n|— >oo V Li J 

Thus, starting from |n| > A, the series Yl, P^) ex P ( — Z±{n)/AKt) in solution (9) can 

|n|=JV v ' 

be majorized by the uniformly convergent series J2 ex P \~ Z±{n) / Ant) . If follows from 

|n|=JV V "'" ' 

here that series in solution (9) are also uniformly convergent. Therefore, we can integrate 
solution (9) term by term and use the majorant estimate 

1 



R(z,z',t) 



(47TKty/ 2 



E 

|ra|=JV 



p(-)(A)exp(-^Up(+)(A)exp^ ^ 



Ant n ' y ' 1 \ Ant 



1 „ (2A- 1)L 

for the remainder of the series which determines Green's function (9). 
It is easy to show that 

Urn *W*o>-*> = 0. 

|n|->oo OZ 
14 



Therefore, the series of derivatives 

oo 

E 



p(-)(Z_(n),A)exp (-^]+ A) exp / Zi(V) 



where 



P« (z±(n)>A) = W.M.a) _ ^M p < ±)(z±(n) A)i 

converges uniformly as well because the series X) pW(Z ± (n), A) exp (— Z±(n)/4«;£) are 

|n|=JV v "' ' 

OO / v 

majorized by the uniformly convergent series — (nL/nt) exp ( — Z±(n)/4/c£J . Thus, so- 

|n|=JV V J 

lution (9) is term-by-term differentiate with respect to z. The possibility of double term- 
by-term differentiation with respect to z and differentiation with respect to t is proved in 
a similar way. 

We represent Green's function (9) in the form 

oo 

G(r,r',t) = £ £ G<?\r,r',t), (A5) 

a=± n=— oo 

where 

eX P Ant 

GW(f, f ', f) = pW(Z ± (n), A) ^ /2 n = 0, ±1, ±2, . . . , (A6) 

substitute expression (A5) into homogeneous equation (1), and carry out term-by-term dif- 
ferentiation of the series. By using recurrence relation (14) for the functions F k (x) , we obtain 
that expression (A5) satisfies the required equations, moreover, each term G^(f, f',t) of 
series (A5) satisfies this equation. By substituting solution (A5) into initial (6) and bound- 
ary (2) and (3) conditions, we make sure that the solution satisfies these conditions. Note 
that, as opposed to the heat conduction equation, boundary conditions (2) and (3) are valid 
not separately for each term G^(f,f',t) of series (A5) but for pairs, namely: 

(§; + A) (GPP, r', t) + G% +1) (f, r>, tj) = 

Thus, relation (9) is the solution of the third boundary-value problem of the heat con- 
duction equation. 
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n = 0,±l,±2,..., 

z=0 

n = 0,±l,±2, .... 

z=L 



[1] H. S. Carslaw and J. C. Jaeger, Conduction of Heat in Solids (Oxford, Clarendon Press, 1959). 

[2] A. V. Lykov and Yu. A. Mikhailov, Theory of Heat and Mass Exchange (Gosenergoizdat, 
Moscow, 1963). 

[3] A. V. Lykov, Theory of Heat Conduction (Vysshaya Shkola, Moscow, 1967). 

[4] A. N. Tikhonov and A. A. Samarskii, Equations of Mathematical Physics (Nauka, Moscow, 
1972). 

[5] N. S. Koshlyakov, E. B. Gliner, and M. M. Smirnov, Partial Differential Equations of Mathe- 
matical Physics (Vysshaya Shkola, Moscow, 1970). 

[6] I. G. Aramanovich and V. I. Levin, Equations of Mathematical Physics (Nauka, Moscow, 1969). 

[7] H. S. Carslaw and J. C. Jaeger, Operational Methods in Applied Mathematics (Oxford, Claren- 
don Press, 1941). 

[8] A. S. Usenko and A. G. Zagorodny, Brownian particle motion in a system with absorbing bound- 
aries, Molecular Physics, 61, No. 5, 1213-1246 (1987). 

[9] M. Abramowitz and I. A. Stegun (eds.), Handbook of Mathematical Functions with Formulas, 
Graphs and Mathematical Tables, (National Bureau of Standards, Appl. Math., Ser. 55, 1964). 

[10] A. G. Butkovskii, Characteristics of Systems with Distributed Parameters (Nauka, Moscow, 
1979). 



16 



